Superdeformed Oblate Superheavy Nuclei? 
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We study stability of superdeformed oblate (SDO) superheavy Z > 120 nuclei predicted by system- 
atic macroscopic-microscopic calculations in 12D deformation space and confirmed by the Hartree- 
Fock calculations with the realistic SLy6 force. We include into consideration high- A' isomers that 
very likely form at the SDO shape. Although half-lives Txi% ^ 10~ 5 s are calclulated or estimated 
for even-even spin zero systems, decay hindrances known for high- if isomers suggest that some SDO 
superheavy nuclei may be detectable by the present experimental technique. 



The question of what is the largest possible atomic 
number Z max of an atomic nucleus is still unsettled. The 
recent experiments on heavy ion fusion in Dubna claim 
Zmax > 118 [1(, with a partial confirmation of hot fu- 
sion cross-sections coming from GSI Q and LBL Berkeley 
Q. Predictions on the stability of superheavy nuclei are 
based either on the Hartree-Fock (HF) studies with some 
effective interaction chosen out of the existing multitude, 
or on the more phenomenological, but also more tested, 
macroscopic-microscopic method. Although these mod- 
els differ quantitatively, they consistently predict prolate 
deformed superheavy nuclei with Z = 100— 112, which 
is confirmed experimentally for nuclei around 254 No 
and spherical or oblate deformed systems with Z > 114 
and N = 174-184, see e.g. 0, ||. In the present letter 
we show that realistic calculations predict superdeformed 
oblate (SDO) nuclei, with characteristic quadrupole de- 
formations —0.4 < /?2o —0.5 (spheroids with the axis 
ratio « 3:2), for Z > 120. By this we confirm and extend 
one of the conflicting conclusions of @ (Fig. 12 there). 
Relying on the calculated energy surfaces, masses and 
cranking mass parameters, we calculate or estimate half- 
lives for selected even-even SDO systems. Then we con- 
sider an idea, advanced e.g. in Q, of extra stable high- if 
shape isomers, also in odd systems, whose existence at 
the SDO shape is very likely. Expected decay hindrances 
point to the possibility that some of these exotic-shaped 
superheavy nuclei, far from the conventionally expected 
"island of stability", live long enough to be detected. 

The Model. Within the macroscopic-microscopic 
method, energy of a deformed nucleus is calculated as 
a sum of two parts: the macroscopic one being a smooth 
function of Z, N and deformation, and the fluctuating 
microscopic one that is based on some phenomenologi- 
cal single-particle (s.p.) potential. A deformed Woods- 
Saxon potential model used here is defined in terms of 
the nuclear surface, as exposed in Q. We admit shapes 
defined by the following equation of the nuclear surface: 

R(0,<p) = c({0})RQ{l + Y,PxoYxo(<B,tp) + 

A>1 

E ^cY^(0,<p)}, (1) 

A>l,/i>0,et;en 

where c({/3}) is the volume-fixing factor. The real-valued 



spherical harmonics Y£ , with even \i > 0, are defined in 

terms of the usual ones as: Yg = (Y\^ + Yx-^)/ V2- 
In other words, we consider shapes with two symmetry 
planes. For the macroscopic part we used the Yukawa 
plus exponential model [9|. All parameters used in the 
present work, determining the s.p. potential, the pairing 
strength and the macroscopic energy, are equal to those 
used previously in the calculations of masses [l(| and 
fission barriers [ll[ of heaviest nuclei. 

Calculations. We used a rich variety of shapes, with 
possible nonaxiality and mass-asymmetry, to reliably de- 
termine energy landscapes of the heaviest nuclei. A de- 
formation set included both traditional quadrupole de- 
formations j3 and 7, where P20 = /3cos7 and /?22c = 
— /3sin7 (for 7 = n x 60°, with n integer, a quadrupole 
shape is axially symmetric), three hexadecapole distor- 
tions /?4o, (342c, /?44c, the higher-rank even axial multipoles 
Peo and f3$o, and the following odd-multipole deforma- 
tions: /330, #320 ^50, /?52c and f3 70 - altogether twelve pa- 
rameters. The range of deformation parameters covered 
a region of shapes up to, and little behind the fission bar- 
rier, where the shape parametrization ([1]) may be hoped 
sufficient. 

Energy landscapes were obtained by a multidimen- 
sional energy minimization on a map of equidistant mesh 
points (/3cos7, flsinj) with respect to 10 other deforma- 
tions. We used a rather large mesh spacing of 0.05 in 
order to make time-consuming calculations feasible. A 
subsequent interpolation served to visualize results. In 
order to check the results we monitored the continuity 
of the resulting 10 deformation parameters with respect 
to /3cos7 and /3sin7, and their stability with respect 
to the choice of their starting values. To assess the lat- 
ter, we repeated the minimization for the whole map for 
selected nuclei by choosing random starting values. We 
have found that the results agreed with the ones obtained 
previously. Additional minimizations have been done to 
further verify the found minima, in particular, the axially 
symmetric minima were reproduced by the minimization 
over the axially symmetric deformations (3xo- 

Equilibria & SDO minima; fission barriers. 
Quadrupole deformations /3 2 o of the g.s. (global) 
minima, calculated for ~ 300 even-even nuclei are shown 
in Fig. Q] In addition to spherical, well- or weakly 
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FIG. 1: Calculated ground state quadrupole deformations P20 
(color online). 



deformed prolate and oblate equilibrium shapes there 
is a region of SDO nuclei for Z > 120, N < 168, of 
particular interest here. SDO global minima occur also 
for large N = 190,192 and Z = 118-122. Although some 
weakly deformed minima have non- axial distortions, 
energies of Z = 120 isotopes plotted vs. P20 for axially 
symmetric shapes in Figj2] fairly illustrate the shape 
competition & coexistence in the Z > 120 region. The 
secondary SDO minima exist there for 168 < N < 172 
and N > 184. They appear also in Z < 118 nuclei. 
Typically, they lie 2 MeV above the g.s. This has 
an effect on the a-decay of the SDO Z = 120 isotopes 
(see below). In the whole Z > 114 region, the deepest 
minima, spherical or oblate, occur for N — 174-184; for 
Z = 124,126 they are predominantly oblate. 

Energy maps in (j3 cos 7, (3 sin 7) plane are necessary to 
appreciate fission barriers, Fig [3] The conspicuous result 
of our calculations is that triaxial saddles occur in all 
studied nuclei. They may lower the axial fission barrier 
by up to 2.5 MeV. This lowering increases with N and 
is larger for bigger Z. The odd-multipole deformations 
do not change the barriers as much, but they lower some 
oblate minima and modify the energy maps around and 
beyond the saddles. 

Crucial for stability is that barriers diminish with N 
decreasing below 174-176 and with Z approaching 126. 
The first feature is common also to the self-consistent HF 
results [l|||, while the second is very distinctive for the 
macroscopic- microscopic model used here (Tlj . Hence, 
the largest barriers of ps 3.4 MeV predicted for SDO nu- 
clei 286 120 and 288 122 are rather small as compared to 
the 5.6 MeV barrier for 296 120 [11]. The barriers for 
N > 190 SDO nuclei are still smaller, so we do not con- 
sider them further. As the a-decay rates increase with 
Z, we concentrate on the SDO nuclei around Z 120 
and TV w 166. 
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FIG. 2: Energy relative to the spherical macroscopic contribu- 
tion, -E(/?2o) — -Emacr (sphere), for the Z = 120 isotopic chain; 
each point results from the minimization over (3\o, A = 3-8 
(color online). 



Other models. To convince ourselves that the SDO 
minima are not a strange twist of the particular model we 
repeated the minimizations for the interesting nuclei by 
using A) the same microscopic model and another version 
of the macroscopic energy, the LSD liquid drop model of 
[lH], B) the selfconsistent HF method with the realistic 
Skyrme SLy6 force [3]. Both calculations support the 
prediction of the global SDO minima; they are even by 
1=3 1 MeV deeper with the LSD variant of the macroscopic 
energy. In the HF calculations, the energy competition 
between prolate, oblate and SDO minima and fission bar- 
riers come out similar as in the macroscopic-microscopic 
study. 

Stability against fission. We checked fission half-lives 
T s f by calculating WKB action with cranking mass pa- 
rameters for selected nuclei. We assumed the zero-point 
energy of 0.5 MeV. To handle fission paths in 12D defor- 
mation space we calculate, instead of the mass param- 
eter tensor, the effective mass parameter along a pre- 
scribed path. Technically, this is done by replacing ana- 
lytic derivatives with respect to deformations by the finite 
differences. 

Two possible classes of fission paths and barriers along 
them may be read from Fig. [3J The barriers along the 
axial saddle (at j3 « 0.3, 7 = 0) are longer and have thin- 
ner peaks. They can compete with the triaxial path only 
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FIG. 3: Energy surface of 288 122, normalized as in Fig[2] 
Crosses mark the saddles (color online). 
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FIG. 4: Mechanism of the a-decay hindrance of the SDO 
286 120; energy normalized as in Fig[2] (color online). 



when there is a deep normal oblate minimum, i.e. for 
A" =166 or 168. Triaxial barriers and the related WKB 
action change smoothly from isotope to isotope. The 
smallest action we found along triaxial, nearly straight 
paths. They give half-lives 1(T 6 s for 286 120 and 10~ 5 s 
for 288 122, with an estimated error of 1 order of magni- 
tude. 

Stability against a-decay. From the calculated masses 
and the improved formula a la Viola-Seaborg [lj| , we ob- 
tain for the g.s.— ?>g.s. transitions Log(T a [s\) = —9.1 for 
286 120, and longer T a for the lighter Z = 120 isotopes. 
These SDO— ^prolate transitions, FigJU must be strongly 



hindered by a very different structure of both configu- 
rations, in particular, the occupation of intruder states 
at SDO shape (see below). If the hindrance would be 
complete, only SDO— >SDO transitions would remain. As 
already mentioned, SDO configurations in the Z = 118 
daughters are excited by 2 MeV (2.5 MeV in HF). This 
leads to a considerable increase in half-life: Log(T a [s\) 
becomes equal to —5.5 for 286 120 and T a are shorter for 
lighter isotopes. Clearly, this result holds for the config- 
uration hindrance factors < 10 -3 6 . 

K-isomers at SDO deformation; odd systems. With 
half-lives T^/ 2 < 10~ 5 s - the present limit for detection of 
synthesized superheavy nuclei - superheavy SDO systems 
might be considered merely as a theoretical curiosity. A 
fascinating possibility for their longer life-times is related 
to K- isomerism, see @, HH ■ Indeed, high- if configura- 
tions at the SDO shape are very likely, see Fig. [5] Due 
to large deformation, the neutron kyj/2 and proton j'15/2 
intruder states with large angular momentum projections 
on the symmetry axis O are close to the Fermi level for 
Z = 120, A~ = 166. Of unique structure and parity, 
they provide identity to high- AT 2(4)-quasiparticle config- 
urations. Candidates for low- lying AT-isomers are the so 
called " optimal" configurations [17j , with singly occupied 
large- O orbitals close to the Fermi level. In 286 120, the 
candidates are the proton (13/2~,7/2 + )10~ and neutron 
(15/2 + ,9/2 _ )12~ configurations. The possible low-lying 
or ground states in odd nuclei are the neutron 15/2+ state 
in 285 120 and the proton 13/2" state in 285 119; the low- 
lying 14~ state could be expected in the odd-odd 284 119. 
Detailed predictions would require energy minimization 
at fixed configuration with blocking. 

In assessing stability of high- AT isomers or odd nuclei 
we rely on estimates and analogies with well established 
experimental facts, as we cannot precisely calculate their 
decay rates. Let us notice that the considered SDO nuclei 
are proton-unstable, but in view of the large Coulomb 
barrier the related life-times may not concern us, at least 
for even- 27 nuclei odd-Z, high- A" states are protected 
by the centrifugal barrier for high-f2 protons. 

Fission hindrance. As well known, T s j for odd and 
odd-odd heavy and superheavy nuclei are by 3-5 orders 
longer than for their even-even neighbours. Similar in- 
crease was found for high-AT isomers, with respect to 
(prolate) shape isomers on which they are built, in even 
240-244Q m p or gjjo superheavy AT-isomers two fac- 
tors combine to increase fission half-life: A) the axial fis- 
sion path is closed by the conservation of the K quantum 
number, B) triaxial barriers increase due to a decrease in 
pairing caused by the blocking of two neutrons or pro- 
tons. Additional hindrance of fission is expected for con- 
figurations involving blocked high-f} intruder states. 

Consider now the effect of the saddle deformation on 
the fission barrier at high spins. The geometrical mo- 
ments of inertia from the HF calculation in 288 120 are: 
J± = 71 b, J|| = 109 b at the SDO shape and j\ = 107 
b at the triaxial barrier. The actual moment of inertia 
at the barrier is reduced by pairing to / ' with f b 
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FIG. 5: Calculated s.p. energies vs projection of the angular 
momentum on the symmetry axis in 286 120 at the SDO g.s.: 
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substantially smaller than 1. Without pairing, J]\ is an 
average moment of inertia of yrast non-collective high- if 
states [13 • As J\\ > f b j\ and pairing at SDO g.s. is 
weaker than at the barrier, there should be no decrease 



with if in fission barrier for SDO if -isomers. 

Alpha-decay hindrance. Although this seems the least 
certain of our arguments, if -isomerism may substantially 
increase a half-lives: the high-if isomer in 270 Ds has 
longer (partial) half-live T a — ms than the g.s., 

T a (g.s.) = 100±4o° ^ El- For SDO nuclei, an additional 
hindrance may result from a difference between the par- 
ent and daughter high- if configuration, or, for the same 
configuration, from its extra excitation in the daughter, 
leading to a smaller Q a . 

Stability against beta-decay. The /3 + decay rates for 
neutron-deficient candidates for the SDO if -isomers can 
be estimated by neglecting the emitted electron energy 
m e c 2 in the decay energy: Qp = (M(A, Z) - M(A, Z - 



1) 



)c 2 . Then, one has \p ~ \M\ 2 G 2 F Q% where |M 



is the transition matrix element and Gf is the Fermi 
constant. Even for a perfect overlap, |Ai| 2 ~ 1, using our 
calculated masses we obtain half-lives Tp — ln2/\p of the 
order of 0.1-1 s for even and odd SDO nuclei, consistent 
with the results by Moller et al. [20]. Since for high- if 
isomers \M\ is reduced, their /3 + decay is even slower. 

Although the production of SDO nuclei is another sub- 
ject, one may notice here that the SDO shape is much 
closer to the sticking point configuration of the prolate 
and spherical heavy ions in the side collision than the 
sphere. 

Summarizing, within both macroscopic-microscopic 
and Skyrme HF methods, one obtains SDO shapes of 
the ground- or low excited states of superheavy Z > 120 
nuclei. Although even-even, spin zero nuclei decay by a 
quick ~ 10~ 5 -10 -6 s fission or a-decay, longer half-lives 
are expected for high- if isomers which very likely exist in 
some even or odd systems. One case of a sizable a-decay 
hindrance could make such a system detectable by the 
present technique. 
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